Realistic Rashba and Dressehaus spin-orbit coupling for neutral atoms 
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We describe a new class of atom-laser coupling schemes which lead to spin-orbit coupled Hamil- 
tonians for ultra-cold neutral atoms. By properly setting the optical phases, a pair of degenerate 
pseudospin (a linear combination of internal atomic) states emerge as the lowest energy eigenstates 
in the spectrum, and are thus immune to collisionally induced decay. These schemes use cyclically 
coupled ground or metastable internal states. We specialize to two situations: a three and four level 
case where the latter adds a controllable Dresselhaus contribution. We describe an implementation 
of the four level scheme for ^'^Rb and analyze its sensitivity to typical laboratory noise sources. 
Lastly, we argue that the Rashba Hamiltonian applies only in the large intensity limit since any 
laser coupling scheme will produce terms non-linear in momentum that decline with intensity. 



Spin-orbit (SO) coupling is essential for realizing topo- 
logical insulators, non-interacting fermionic systems with 
topological order [H and yet in other contexts it 
leads to parasitic effects such as reduced spin coher- 
ence times [3]. As with the progression from the single- 
particle integer quantum Hall effect (IQHE) to the inter- 
action driven fractional quantum Hall effects (FQHEs), 
the next important step is realizing the strongly interact- 
ing cousins to the topological insulators, of which topo- 
logical superconductors are a first example 4, 5J. Since 
ultracold atoms lack intrinsic SO coupling, numerous 
techniques for generating SO coupling (generally equiva- 
lent to nonabelian gauge potentials [5]) with optical [3- 
[TT| and now rf [12 fields have been suggested, one of 
which was recently implemented |13j . 

Current proposals for realizing SO coupling suffer from 
two primary limitations. First, the pair of dressed spin 
states comprising the effective spin- 1/2 system are not 
the two lowest energy states, so coUisional deexcita- 
tion [14] can rapidly transfer population into the ground 
state [15]. Second, the required tripod coupling scheme 
is difficult to directly realize in alkali atoms [TB]. In this 
Rapid Communication, we introduce a class of laser cou- 
pling techniques that overcome these difficulties, and we 
explore the departure of such models from the ideal case 
(suitable only at infinite laser intensity). 

In condensed matter systems, SO coupling links the 
linear or crystal (not orbital) momentum fi,k to the spin 
of an electron, for example. For systems confined to two 
dimensions (2D), terms linear in momentum can be rep- 
resented a,s cL sum of R,8-sliba, o. {&x^y — ^ykx ) and Dressel- 
haus /3 {a^ky + CTyfca;) SO couplings [I71[TH|, where d'x,y,z 
are the Pauli matrices. Proposals for creating SO cou- 
pling with neutral atoms use lasers to link states of dif- 
ferent momentum and spin. Because these lasers impinge 
from discrete directions, the system lacks the continuous 
rotational symmetry of the pure Rashba Hamiltonian an- 
ticipated in earlier works [TUl [HI [H] . We show how a 



perturbative treatment restores the system's N-fold ro- 
tational symmetry and demonstrate that these pure cou- 
plings are only exact in limit of infinite laser-atom cou- 
pling strength. 

We consider N ground or metastable atomic 
"spin" states {|1) , |2) , • • • , |iV)} cyclicly coupled to- 
gether with complex valued matrix elements f2j+ij- — 
—fl exp [i (kj -x -|- 7j )] representing plane waves that link 



(a) Coupling scheme (b) Geometry 




(c) Uncoupled dispersion (d) Coupled dispersion 




FIG. 1. Four level scheme, (a) Effective coupling between 
four internal ground or metastable atomic levels, (b) Spatial 
orientation of coupling fields, (c) Uncoupled eigenenergies for 
fl = e = 0. The four free parabolae are displaced by Kj 
in the Qx-Qy plane, (d) Dispersion of four dressed states for 
n = SEr, e = 0, and 7 = 7r/4, showing the expected Dirac 
points, one for each pair of dressed bands. 
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consecutive states \j) to + Here, describes the 
optical coupling strength; Kkj and jj are the respective 
discreet momentum and phase acquired in the j -> j + 1 
transition. Throughout this manuscript we apply "peri- 
odic boundary conditions" \N + 1) = |1) for spin states. 
See, for example, the 4-state topology in Figure [l^. 

Including the motional degrees of freedom, the many- 
body Hamiltonian 



a. Three level: Af = 3, 7 = vr/S b. Four level: W = 4, 7 = 7r/4 
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describes a system of atoms with mass ra in 2D absent 
the ubiquitous confining potential. Here, {(/)^(k)} is the 
spinor field operator describing the creation of an atom 
with momentum Sk in internal state for even N ^ we 
introduce e, describing a detuning of alternating sign. 

In what follows, we require that ^ k.^ = 0, so that 
no momentum is transferred to an atom during a |1) — 
|iV) — )■ |1) transition. We define the momenta- 



FIG. 2. Eigenenergies Et of the coupling term in Eq. S 
showing the cosinusoidal energies evaluated at integer I. (a) 
Three level case with 7 = 7r/3 and (b) Four level case with 
7 = 7r/4. For both cases we identify pseudospin states |t) and 
14,) with the lowest energy pair of states. 

for p G Z. The manufactured degeneracy between states 
at £ = and 1 for 7 = 7r/iV is illustrated by Fig. [2] for 
iV = 3 and 4. 

When the displacement vectors Kj reside on the ver- 
tices of a regular polygon, Kj — — ki^ sin {27Tj /N) + 
ki, cos {2Trj/N) Sy, the full Hamiltonian matrix is 



exchange with differences kj = Kj_|_i — and require jj^ (^q-j _ i §^ 



{Kj} to have zero average. Moreover, 7 may replace the 
phase 7j of each state vector where 7 = ^^li/N 
PP] , without loss of generality. 



With the substitution (^sHq) 



K, 



the Hamiltonian [Eq. ([T])] separates into an integral 
/Ej,j' <^](q)^^j'(q)<^j'(q)rf'q/(27r)2 over iVx^ blocks 



- , , /i^ |q + K,-l / e c 
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each labeled by a quasi-momentum /iq. The first (ki- 
netic energy) term in Eq. ([2| describes the 2D displaced 
parabolae depicted in Fig. In analogy to band- 

structure, the last (coupling) term in Eq. ^ has the form 
of an N site ID periodic tight binding Hamiltonian with 
a "magnetic flux" A^7 and a "hopping" matrix element 
il/2, where internal atomic states play the role of lattice 
sites. We diagonalize the coupling term, (a zero-order 
approximation suitable when 51 is much larger than all 
other parameters), by transforming into the basis conju- 
gate to the spin-index j with field operators 



+ [{iqx + Qy) + h-C] + -5f-JV/2,f', (3) 

where momenta and energies are expressed in recoil units, 
fci and E'l = h?k\/2m respectively. Assuming 7 « 7r/iV, 
we focus on the manifold of two nearly degenerate states 
with £ = and i = 1, yielding the pseudospins jj,) and 
It) depicted in Fig. [2) 

In what follows, we derive an effective 2x2 Hamil- 
tonian H sa + (2) + 7/(3) foj. tiiis manifold up to 
third order in powers of In the subspace spanned 

by the lowest energy pseudospin-pair we obtain (up to a 
constant) a zero-order Hamiltonian of the Rashba form 



= |q|' i + {a.,qy - dyq,) 
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The diagonalization provides the eigenenergies Ei = 
— rj cos(27r£/iV — 7) of the coupling Hamiltonian, where 
£ G {0, - ■ ■ , iV — 1} is analogous to the usual crystal mo- 
mentum. The ground state can be made two-fold degen- 
erate by tuning 7 to "magic" phases 7 = 2tt{p + 1/2) /N 



with a Zeeman field lS.z — Ei — Eq ^ —2517' sin (vr/A^) 
generated by slight changes 7' — 7 — tt/A from the magic 
phase. At finite coupling SI, we adiabatically eliminate 
the excited states order-by-order in perturbation theory 
giving effective terms TJ*^"-' in the ground manifold Hamil- 
tonian. 

Since Eq. (|3]) is cyclic for e = 0, we expect an energy 
shift at order n = 2 in perturbation theory (effectively 
a Stark shift), and pseudospin-changing terms at order 
n = A — 1. These terms serve to restore the A^-fold rota- 
tional symmetry absent from Eq. Q and in the analogous 
expressions of earlier proposals [HHHlIin]- To understand 
the departure from the Rashba Hamiltonian, we will first 
consider the simpler e = case. 

For the A^ = 3 case the second order effective Hamil- 
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tonian 



(a) 4-Level dispersion at SEi coupling 
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|q|'i- 



restores the expected 3-fold symnietry. For iV = 4, the 
second order term acts as a state- independent Stark shift, 
but the third order correction 
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[^x{ql - iqyql) ~ ^y{ql - Sg^g^)] 



restores the 4-fold rotational symmetry. This term is 
reminiscent of the cubic Dresselhaus SO coupling present 
in GaAs 2D electron systems [21 HI] ■ The Rashba Hamil- 
tonian's [Eq. Q] ground state energy is minimized on 
the ring where |q| — 1/2; the perturbative terms mod- 
ulate both the momenta and energy where the minima 
occur [22]. Figure[3^ shows this modulated energy for the 
4- level case at f2 = 3i?i,; Fig. [3]d plots the peak-to-peak 
amplitude of the energy modulations. The perturbative 
analysis given by the dashed lines rapidly converges to 
the numerical solution to the full Hamiltonian given by 
the solid curves. Since the magnitude of the modulations 
scale like 17^"^, the 3-level case requires an impractically 
large f2 > IOQEl to reduce the corrections to the Rashba 
Hamiltonian below 10~^El, while the 4- level case re- 
quires just r2 = IOEl. Interestingly, the familiar tripod 
scheme [SI |H1 HOI Ell ES] reduces to our = 3 ring model 
when far detuned from the excited state, with 7 = for 
red detuning and 7 = 7r/3 for blue detuning. By contrast, 
the corrections to the Rashba model converge as in 
the standard resonant tripod (a four level system). 

The alternating detuning e featured in Eq. ^ leads 
to a super-lattice in the above mentioned band-structure 
analogy, and its contribution can be included exactly. 
However, for a more painless description, we take e/fi <C 
1, which for iV = 4 adds a tunable Dresselhaus term 
Hb ~ (^ELe/y/2ilkL^ {a^qy + o'yqx) at second order in 
fl^^ (in original units). Thus, our scheme produces 
both Rashba and Dresselhaus couplings with strengths 
a — E^/kL and (3 — Ei^tj\f2^kL along with a z-aligned 
Zeeman field Az « — v^i^T"'. The laser configuration 
specifies a; the alternating laser detunings set /3; and Az 
is controlled by the average phase 7. Together this allows 
for nearly complete control of the SO Hamiltonian. 

Implementation In the following, we specialize to a 4- 
state configuration within the F = 1 and F = 2 hyperfine 
manifolds of ^^Rb's ^Si/2 ground electronic state (illus- 
trated in Fig|4^). Each pair of states is Raman coupled 
by two lasers tuned between the ^P3/2 and atomic 
excited manifolds {El « /i x 2500 Hz). Although there 
are nominally eight transitions coupling four states we 
reuse three coupling lasers. The symbols ai, 02, ca, tti 
and 7r2 each implicitly denote a laser frequency (that may 
be inferred from the level diagram in Fig|4^, the energy of 
each state | j) and the choice of 2-photon detuning for one 
set of Raman lasers) and a direction given by Figl4|D. The 
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(b) Scaling 
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FIG. 3. a) Energy of the lowest eigenstate for N = 4, 
Q, = S-El, e — 0, and 7 = 7r/4 plotted in the qj:-qy plane 
showing the four-fold rotational symmetry. The graph is col- 
ored according to the energy, and white curves mark contours 
of equal energy. The yellow dashed line depicts azimuthal 
modulation of radius in momentum of the energy-minimum, 
b) Peak-to-peak magnitude of the azimuthal energy modu- 
lations -Emod plotted as a function of coupling fi, showing a 
rapid suppression for A*' = 4 as compared to A^ = 3. The 
solid red curves, the result of exactly solving the full Hamil- 
tonian, quickly converge to the perturbation result (black 
dashed lines). 

first three symbols denote left circular polarization and 



their relative phases define A7 — 20^ 



Lasers 



TTi and 7r2 are linearly polarized along the magnetic field 
and produce canceling phases in the Hamiltonian and 
thus do not contribute to the phase-sum. Additionally, 
state labels will match those of the paper given order- 
ing of frequency ly^^^ > v„.^ > v^,^ > > A bias 
field along ij and an ac Zeeman shift provided by a mi- 
crowave field along x produce a composite Zeeman shift 
that spectroscopically isolate JJ! the hyperfine changing 
transitions by at least lOfi. In addition, |2) and |4) shift 
together with applied magnetic field, allowing straight- 
forward control of the detuning e, thus contributing Dres- 
selhaus coupling. 

In practice, electronic equipment will introduce un- 
wanted time-varying terms into the Hamiltonian. How- 
ever, the effects of magnetic field detuning Ae, laser in- 
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FIG. 4. 4-level ring coupling scheme in *^Rb involving 
hyperfine states ttif) Raman-coupled by a total of five 
lasers marked cri, (T2, (Ts, tti and 7r2. (a) Level diagram; 
Each red curve represents a two-photon Raman transition 
with polarizations as marked by the field-symbols. The ring 
states are mapped to physical states according to |1) = |2, 0), 
i2) = 11,-1), |3) = |1,0), and |4) = |2,1). A 6.8 GHz mi- 
crowave field with coupling strength fi^ ac Zeeman shifts the 
internal states |1,0) and 12,0). (b) Schematic layout: a bias 
magnetic field Bo ~ 0.2 mT lies along y. Lasers tti and tv^ are 
TT polarized. Lasers ai, a2 and (J3 are left circularly polarized 
and their relative phases define A^7 = 200-2 ^ 4'cri — (^t^- 



The proposed coupling scheme provides a robust plat- 
form for generating SO coupling for neutral atoms. Be- 
cause the two spin states are the lowest energy dressed 
states, atom-atom interactions cannot induce coUisional 
decay [HI [TS]. In addition, this technique can require 
considerably less laser intensity than prior far-detuned 
schemes to reach nearly pure SO coupled Hamiltonians, 
greatly reducing spontaneous emission. Lastly our spe- 
cific implementation uses only the Smp = 0, ±1 Raman 
transitions allowed at large atomic detuning in the alkali 
atoms p[ ITT]. 
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tensity shifts Afi and phase shifts A7 can have negligi- 
ble contributions in typical laboratory settings. In the 
scheme shown in Fig. |4j time varying magnetic fields in- 
duce identical linear Zeeman shifts of |2) and |4) that 
cancel at lowest order in perturbation theory (while |1) 
and 1 3) are first order insensitive to magnetic fields). 
These Zeeman shifts give an additional Dresselhaus con- 
tribution H/^D cx Ae {axqy + dyqx)', near the minima of 
the Rashba Hamiltonian, the energy shift from the re- 
sulting Dresselhaus term would be 0.0005£'l for our 
Ae « X 50 Hz detuning noise amplitude [T3] . 

Shifts (resulting from laser intensity noise) in a sin- 
gle fij.j+i matrix element add coupling terms 7?An ~ 
— V^Afij .,+i(itj +Kj-^x) ■ o'/Sfci, where & is the vector 
of Pauli matrices. For lasers stabilized at the 0.1% level, 
and with — 10i?L, the unwanted coupling terms have 
magnitude 0.002^;^. 

Phase shifts A7 directly modulate the effective Zee- 
man fields -ffA0 ~ —^/2^l/^'y^z/'2' which opens a gap at 
the Dirac point. For ^^Rb, we require the resulting gap 
energy to be less than a typical T = 10 nK temperature. 
For Q, — 10i?L, the RMS phase noise must have an am- 
plitude A7 < 0.04 rad. Even for independent lasers, this 
level of phase stability is routine [23] . 

Discussion Inspection of Eq. ([s]) shows that Dirac 
points are present for all and e provided 7 = 7r/iV, 
even when perturbative corrections are important. As a 
result, properties of fermion systems which depend only 
on the topology of the dressed-state dispersion may be 
insensitive to small corrections to Eq. (|3]). In contrast, 
bosons generally condense at the energy minima. Thus, 
for insufficient laser coupling, local minima may spoil 
correlation physics potentially arising from the Rashba 
Hamiltonian's degenerate ring of minima [191 124] . 
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